Entangled spin clusters : some special features 
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Abstract 

In this paper, we study three specific aspects of entanglement in small spin clus- 
ters. We first study the effect of inhomogeneous exchange coupling strengths on the 
entanglement properties of the S=^ antiferromagnetic linear chain tetramer com- 
pound NaCuAsOi. The entanglement gap temperature, Te, is found to have a 
non-monotonic dependence on the value of a, the exchange coupling inhomogeneity 
parameter. We next determine the variation of Te as a function of S for a spin dimer, 
a trimer and a tetrahedron. The temperature Tg is found to increase as a function of 
S but the scaled entanglement gap temperature Ie goes to zero as S becomes large. 
Lastly, we study a spin-1 dimer compound to illustrate the quantum complementarity 
relation. We show that in the experimentally realizable parameter region, magneti- 
zation and entanglement plateaus appear simultaneously at low temperatures as a 
function of the magnetic field. Also, the sharp increase in one quantity as a function 
of the magnetic field is accompanied by a sharp decrease in the other so that the 
quantum complementarity relation is not violated. 



I. INTRODUCTION 

Entanglement is a key feature of quantum mechanical systems and gives rise to non-local 
correlations over and above those expected from classical considerations [l]. It can be 
of different types : bipartite, multipartite, zero-temperature, finite-temperature etc. for 
which suitable measures are available in certain cases OEimiHlEllZllHl- In the past 
few years, quantum spin systems have been extensively studied to gain knowledge on the 
different aspects of entanglement. The spins in such systems interact via the exchange 
interaction and also with an external field, if any. Several studies show that the amount 
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of entanglement can be changed by varying the temperature T and/or the magnitude of 
the external field [Hllllini- In the case of entangled thermal states, one can define a critical 
temperature below which entanglement is present in the system and above which entangle- 
ment vanishes, i.e., the system becomes separable. Detection of entanglement can be made 
with the help of an entanglement witness (EW) which is an observable the expectation 
value of which is positive in separable and negative in entangled states [IHl IHl Il2] . Ther- 
modynamic observables like internal energy, magnetization and susceptibility have been 
proposed as EWs O [13 IH] . There is now experimental evidence that entanglement can 
affect the macroscopic properties of solids like specific heat and magnetic susceptibility 
Recently, it has been shown that for separable states, the sum of magnetic susceptibilities 
in the three orthogonal directions x, y, and z obeys the inequality 



X — Xx ~l" Xy ~l~ Xz 



> 



NS 



(1) 



where is the total number of spins in the system, S the magnitude of the spin, ks the 
Boltzmann constant and T the temperature If the magnetization operator = 

J2j commutes with the Hamiltonian H of the system, i.e., [H,Ma] = 0, the magnetic 
susceptibility Xa =x, y, z) can be written as 



Xc 



(2) 
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Thermodynamic properties in general relate to macroscopic systems and the thermal state 
of such a system is entangled if x is < (Eq. (1)). Using the susceptibility inequality 
as an EW, one can detect entanglement from the experimental data without requiring a 
knowledge of the Hamiltonian of the system. 

For a multipartite Hamiltonian H, one can define the entanglement gap as 



E - E 

-'-'Sep 
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where Esep is the minimum separable energy and Eg the ground state energy of the Hamil- 
tonian [12]. For a spin Hamiltonian, Egep is the ground state energy of the equivalent 
classical Hamiltonian [TT]. If a system has entanglement gap G^; > 0, then one can define 
the entanglement gap temperature, T^, as the temperature at which the thermal (internal) 
energy U{Te) = Egep- For temperature T < Te, the thermal state of the system is bound 
to be entangled. Recently, a quantum complementarity relation has been proposed [14] 
between the thermodynamic observables, magnetization and magnetic susceptibility. This 
is given by 



NS N-^S^ 



< 1 



(5) 
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where (US^ = {M^f + {Myf + {M^f. Define the quantities 



p = \ L = 1- ^ (6) 

The quantity P, which depends upon the magnetization, describes the local properties of 
individual spins whereas Q, which involves the susceptibility, is representative of quantum 
spin-spin correlations. From Eq. (1), a nonzero positive value of Q implies the presence 
of entanglement in the system, i.e., non-local correlations. The complementarity relation 
shows that the non-local properties are enhanced at the expense of the local properties in 
order that P + Q is < 1. 

The EWs based on the internal energy and susceptibility have been used to study the en- 
tanglement properties of the spin-| antiferromagnetic (AFM) compounds Cu{N03)2, 2.51)2 
0(CA^)(system of weakly coupled spin dimers) (NHEt).^ [vi^V/ AsgO^o {H2O)] .H2O 
(system of weakly coupled spin tetramers) [in] and the nanotubular system Na2V307 (con- 
sists of weakly-coupled nine-spin rings) |T2|. The weak coupling between the spin clusters 
allows each system to be treated as consisting of effectively independent clusters. Since the 
clusters contain a few spins, the theoretical calculation of entanglement-related quantities 
becomes possible. A number of molecular magnets are known which are well-described in 
terms of small spin clusters such as dimers, trimers, tetramers, tetrahedra etc [IH]- For 
non-bipartite clusters with "all-to-all" spin couplings (trimers, tetrahedra), the EWs based 
on the internal energy and susceptibility give the same estimate of the temperature above 
which entanglement vanishes [H]. For bipartite clusters (a tetramer describing a square 
plaquette of spins with only nearest-neighbour (NN) exchange couplings provides an ex- 
ample), the EW based on the internal energy can detect only the bipartite entanglement 
between two qubits |12|. The spin clusters considered so far are described by Hamiltonians 
with homogeneous exchange interaction strengths. In this paper, we consider the S = ^ 
AFM linear tetramer compound NaCuAsO^ in which the linear tetramer consisting 
of four spins is described by the Heisenberg Hamiltonian 

Hlt = JSl.S^ + aJS^.sl + JSl.sl (7) 

We study the entanglement properties of this compound using both the internal energy 
and the susceptibility as EWs. We next determine the entanglement gap temperature Te 
of small spin clusters as a function of the magnitude 5" of spins. Lastly, we determine the 
quantities P and Q (Eq. (6)) appearing in the complementarity relation (Eq. (5)) for a 
spin-1 dimer compound [Ni2 {Medpt)2{^—ox){H20)2\{ClOi)2-'2.H20 [213] and show that the 
sharp changes in the magnetization and the formation of plateaus at low temperatures are 
accompanied by sharp changes and plateaus in the amount of entanglement. Magnetization 
plateaus have been observed experimentally in the spin-1 dimer compound. This compound 
thus provides a concrete example of a system in which the amount of entanglement can 
change steeply as a function of the magnetic field or does not change over a range of field 
values. 
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FIG. 2. Plots of Tq^ and Tq^ versus a. Tq is the critical entanglement temperature for 
the pair of spins at sites k and I. 
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11. LINEAR CHAIN TETRAMER 



The 5* = I AFM compound NaCuAsO^^ has a linear chain tetrameric structure described 
by the Hamiltonian, -^lt, in Eq. (7) with a ~ 0.4. The term "linear" refers to the pattern 
of exchange couplings and not to the spatial structure of the tetramer [121. The total spin 
5"*°* of the tetramer has the values 2, 1 and 0. There are five S"*"* = 2 states, nine 5**°* = 1 
states and two 5**°* = states. The different eigenvalues and eigenvectors are displayed in 
Appendix A. We now discuss the finite-temperature entanglement properties of the linear 
chain tetramer. The thermal density matrix p (T) is given by 

pm = 4 EEe-'^''M^.-)(^.™l (8) 

^ m 

The first summation is over all the independent energy eigenstates and the second sum- 
mation includes terms corresponding to the (2S'*°* + 1) degenerate eigenstates with the 
eigenvalue E^. Z denotes the partition function 

Z = Y. (25*°* + l) e-^^^ (9) 

A measure of entanglement between the spins at sites k and / is given by the concurrence 
Cki- This is calculated from the reduced thermal density matrix pki (T) using standard 
procedure |2lin|. Since the eigenvalues and eigenvectors of the linear chain tetramer are 
known, the calculation of Cm is straightforward. Figure 1 shows the variation of Cu as 
a function of temperature for a = 0.4 and ^ = 92.7 K, the parameter values relevant 
for NaCuAsOi [IHllin]. The concurrence C23 is zero for these parameter values. One can 
further define a critical temperature T^^ above which the entanglement between the spins at 
the sites k and I disappears. Figure 2 shows a plot of and T"^ versus a for ^ = 92.7 K 
. We next calculate the entanglement gap temperature Te at which the internal energy 

where Egep is the minimum separable energy. Figure 3 shows a plot of Te versus a for ^ = 
92.1 K (Curve a). The operator H - E^ep is an EW since Tr [p {H - Esep)] = U (T) - E^ep 
is < (> 0) when the thermal state is entangled (separable). If the Hamiltonian H of 
the system contains only local interactions such that H = J2<ij> Hij and the underlying 
lattice is bipartite, then 

H Eg^p ^ ^ {.-^ij ^sep,ij) (-i-^) 
<ij> 

where esep,ij is the minimum separable energy associated with the interaction between the 
spins located at the sites i and j. In the case of a translationally invariant Hamiltonian, 
Hij and esep,ij are the same for each interacting spin pair. Each term in the sum on the 
RHS of Eq. (11) can be considered as a bipartite EW. Thus the expectation value of 
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H — Esep is negative only if the two spins in the interacting spin pairs are entangled. In 
the case of the linear chain tetramer, the Hamiltonian Hlt is not translationally invariant. 
This is true even in the limit a = 1. Since, the Hij's and esep,ij's are no longer the same 
for each interaction bond, T^? = T^" ^ T^?. The expectation value of H — Egep now 
depends on the relative magnitudes and signs of the two types of terms on the RHS of Eq. 
(11). In contrast, consider a closed chain of four spins in which the NN spins interact with 
the same exchange interaction strength. In this case, because of translational invariance, 
T}? = Tq^ = T^^ = = Tc and the entanglement gap temperature Te is equal to Tc, 
the critical temperature beyond which the entanglement between two NN spins vanishes 
(concurrence is zero). In the case of the linear chain tetramer, a similar interpretation 
cannot be given. 

We now use the magnetic susceptibility as an EW to determine the critical temperature 
beyond which the thermal state of the linear chain tetramer is separable. We consider 
the case of zero-field susceptibility. In the absence of a magnetic field, (Mq.) = (a = 
x,y,z). Also, due to the spin isotropy of the Hamiltonian, Hlt, (5'*°* is a good quantum 
number), Xx = Xy = Xz = X- The susceptibility x can be written as 



X = 7^ E (S'^*"* + l) {S'"' + l) S'^'e-"''^ (12) 



The susceptibility inequality for separable states (Eq. (1)) becomes 

NS , , 

The critical temperature Tq is given by the intersection point of the two curves : x versus 
T plot from Eq. (12) and x versus T plot from the equality in Eq. (13) [HI HE]- For 
a = 0.4 and ^ = 92.7 one obtains the estimate = 90.88 which is really the 
lower bound of the critical temperature above which entanglement vanishes. Figure 3 also 
shows the variation of as a function of a (Curve b). Since the linear chain tetramer is 
associated with a bipartite graph, is > Tg;, the entanglement gap temperature. 

III. GENERAL SPIN S AND Te 

We have so far considered the case 5' = |. We now consider dimers, trimers and tetrahedra 
of spins of magnitude S. The Hamiltonians describing the small clusters are 

Hdimer = J S1.S2 (14) 
Htrimer J 
H tetrahedron J 
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FIG. 3. Plots o^Te (Curve a) and (Curve b) for the linear chain tetramer as a function 
of a (^ = 92.7 K). 



A trimer and a tetrahedron are defined on a non-bipartite graph. The AFM cluster Hamil- 
tonian in each case contains "all-to-all" spin couplings and is frustrated as there is no 
separable state that simultaneously minimizes the energy of each interacting spin pair. 
Thus esep,ij in Eq. (11) is no longer the minimum separable energy of an interacting spin 
pair, it has a magnitude greater than that of the latter quantity. The minimum separable 
energy for the whole Hamiltonian is Egep = J2<ij> ^sep,ij = -^tot e^ep, where Ntot is the total 
number of interacting spin pairs. Since e^ep is greater than the minimum separable energy 
for an interacting spin pair, the EW, H — Eg^p, can detect entanglement even if the en- 
tanglement between the spins in the interacting spin pair vanishes, i.e., the corresponding 
reduced density matrix becomes separable [12]. In this case, the entanglement gap tem- 
perature Te is Tc, the critical temperature beyond which the NN concurrence is zero. As 
shown in |l^, Te = Tq ior S = ^ non-bipartite clusters like the trimer and the tetrahe- 
dron described by the Heisenberg Hamiltonian with homogenous exchange couplings. This 
result holds true for general S in the case of spin clusters with "all-to-all" homogeneous 
Heisenberg spin couplings. We thus use only the internal energy-based EW to determine 
how the critical entanglement temperature varies as a function of S in the cases of the spin 
dimer, trimer and the tetrahedron. 

For Hamiltonians with "all-to-all" spin couplings, the energy eigenvalues of all the eigen- 
states can be determined quite easily from a simple formula. The Hamiltonian can be 
written as 




(17) 



where 



tot 



J2iLi The eigenvalue Estot for a state with total spin 5'*°* is 
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Eqtot = — 

* 2 



where S is the magnitude of a spin. The possible values of 5"*°* are NS, NS — 1, etc. 

The lowest value is zero for N even and | for N odd. Under the vector addition of angular 
momenta, a particular 5"*°* value can be achieved in more than one way, i.e., has some 
multiplicity. Let -Pftot^v be the multiplicity, i.e., the number of possible states with total 
spin angular momentum 5**°* when spins, each of magnitude S, are combined. As shown 
by Mikhailov -PftotAr is given by 

^ E (-1)' ( : ) ( + ^' - " - ' ) (IS) 

CTTl \ 
are the binomial coefficients. The summation index k satisfies two conditions 
n J 

: (i) k > and (ii) the upper numbers in the binomial coefficients cannot be less than 
the lower numbers. Thus, < k < [ ^^^s+i ^ J '^here [b] denotes the integer part of b. 
The minimum separable energy of a spin cluster is equal to the ground state energy of 
the equivalent classical Hamiltonian. In the classical ground state, 5"*°* = and each 
< Sf >= S"^. Thus, the minimum separable energy, Eg^p, for the dimer, trimer and 
tetrahedron is given by Eggp = — 5^ (dimer), — (§5'^) (trimer) and —25^ (tetrahedron). 
The entanglement gap temperature Te can be calculated by using the relation in Eq. (10). 
Figure 4 shows the variation of Te with S for dimers (star), trimers (solid square) and 
tetrahedra (solid diamond). The entanglement gap temperature, Te, is found to increase 
with S in each case. According to conventional notion, spins behave as classical objects in 
the limit of large S. The commutation bracket of spin operators, with each operator scaled 
by the total spin S, tends to zero as S' — > cxd. One would thus expect the entanglement gap 
temperature Te to decrease rather than increase as the magnitude of S is raised. Some 
earlier studies have reported findings similar to ours. Hao and Zhu |22] have studied the 
AFM Heisenberg chain with spins of magnitude S. For a two-sited chain, i.e., a dimer, 
they find that the entanglement gap temperature Te increases almost linearly with S. 
For S = 1, they have shown that Te decreases as the length of the chain is increased. 
Wiesniak et al. have determined the critical entanglement temperature, T^, based on 
the susceptibility as an EW, and find the result that Tq = 1.6 J for the S = ^ Heisenberg 
chain and = 2 J for a chain of spins 1. As pointed out by Bowling et al. [12], it is 
sensible to define a scaled temperature 

ksT , , 

t = (20) 

for a meaningful comparison of Hamiltonians with different total energy ranges, Etot {Etot 
is the difference between the highest and the lowest energy eigenvalues). The scaled en- 
tanglement gap temperature can be defined as Ie = The inset of Fig. 4 shows 
the variation oitE with S for dimers (star), trimers (solid triangle) and tetrahedra (solid 
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FIG. 4. Variation of the entanglement gap temperature, Tg, with S for dimers (star), 
trimers (sohd triangle) and tetrahedra (solid square). The inset shows the variation of the 
scaled entanglement gap temperature with S. 

square). One finds that the scaled entanglement gap temperature decreases as S increases. 
The result can be interpreted in the following way. As S increases, the fraction of the 
total energy range of the spin system which corresponds to entangled states decreases and 
tends to a limiting value as S" — oo. Classical behaviour presumably emerges when the 
entangled states have a negligible contribution to the total energy range. 



IV. SPIN-1 DIMER : QUANTUM COMPLEMENTAR- 
ITY 



Bose and Chattopadhyay |23 have considered some toy spin models and shown that first 
order quantum phase transitions, occurring at special values of the external magnetic field, 
are accompanied by magnetization and entanglement jumps. Upward jumps in the magne- 
tization give rise to downward jumps in the amount of entanglement. Also, magnetization 
and entanglement plateaus coexist in the same range of magnetic fields. Later studies 
established the general validity of these results [211123 EE] • In this section, we show that 
the quantum complementarity relation (Eq. (5)) provides a natural explanation for the 
correlated changes in the amounts of magnetization and entanglement as a function of the 
magnetic field. We illustrate this in the case of a spin-1 dimer compound [Ni2 {Medpt)2{iJ'— 
ox){H20)2\{Cl0^2-'^H20 [Medpt = methyl — bis{3 — aminopropyl)amine) which exhibits 
magnetization plateaus at sufficiently low temperatures [2il|. The Hamiltonian describing 
the spin-1 dimer is 



H, = J {SfS^ + SfSl) + 5 J (S^SI) + d {Slf + {Slf + B {SI + ^2^) (21) 
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FIG. 5. Plots of P, Q and P + Q as a function of j in the case of a spin-1 dimer 
(S=l,d = Q and (3J = 3.) 

where 6 is the exchange anisotropy parameter, d labels the axial zero-field splitting pa- 
rameter and B the strength of the external magnetic field. The negative (positive) sign of 
the parameter d corresponds to an easy-axis (easy-plane) single ion anisotropy. If the spin 
system is entangled, a sharp increase in the magnetization (obtained at low temperatures) 
is accompanied by a sharp decrease in the amount of entanglement so that the complemen- 
tarity relation is not violated. At T = 0, the sharp changes become the 'jumps' associated 
with first order quantum phase transitions. As T increases, the changes occur more gradu- 
ally as a function of the magnetic field. In an entangled system, the magnetization plateaus 
are accompanied by entanglement plateaus and the complementarity relation continues to 
be valid. To illustrate this, we first calculate the eigenvalues and the eigenvectors of the 
dimer Hamiltonian Ha (Eq. (37)). These are displayed in Appendix B. The magnetization 
M (only the z-component is non-zero) and Xz, the z-component of the susceptibility are 
derived from 

1 dZ dM 

The susceptibility components Xx and Xy are determined from Eq. (3) with (S^) , {Sf) , {S2) 
and {S2) — since the magnetic field is in the z-direction. One can now calculate the terms 
P and Q (Eq. (6)) appearing in the quantum complementarity relation given by Eq. (5). 

Figure 5 shows the plots of P, Q and P -|- Q as a function of -j for (5 = 1, = 
and PJ = 3. Figure 6 shows the appearance of plateaus as the temperature is lowered 
(/3J = 20). Note that a sharp increase in P is accompanied by a sharp decrease in Q. 
Plateaus in P and Q occur in the same range of magnetic fields. At sufficiently low 
temperatures, the two-step plateau structure is still obtained for non-zero values of d. The 
intermediate plateau has a lesser width for negative values of d and disappears at d — —1. 
At this point, Q is < throughout the range of -j values so that the spin system is not 
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FIG. 6. Plots of P, Q and P + Q as a function of -j in the case of a spin-1 dimer 
(S = l^d = and pJ = 20.) 

entangled. If d is changed from d = —1 to d = +1, the two-step structure in both P and Q 
is recovered and the amount of entanglement is no longer zero. The easy-plane single ion 
anisotropy {d > 0) is found to be favourable towards plateau formation in both P and Q. 
The changes in P and Q are correlated so that the complementarity relation P + Q < 1 
is always valid. A two-step magnetization curve has been experimentally observed in 
the spin-1 nickel compound mentioned earlier [201 . Theoretical calculations, based on a 
description of the compound as a collection of independent spin-1 dimers, give a good fit to 
the experimental data on the magnetization and susceptibility. The exchange anisotropy 
parameter 6 has been taken as 1 and the single ion anisotropy is of the easy-axis type 
{d < 0). Magnetization experiments have been carried out for the external magnetic field 
parallel to the z and x directions. In both the cases, a two-plateau structure has been seen 
in the magnetization versus field curves. In the latter case, the plateau structure is found to 
be more prominent. Our theoretical calculations suggest that the magnetization plateaus 
exhibited by the spin-1 nickel compound are accompanied by entanglement plateaus. 

V. SUMMARY AND DISCUSSION 

In this paper, we study some special features of entangled small spin clusters. We first 
consider the S = ^ AFM linear chain tetramer compound, NaCuAs04, described by 
the Heisenberg exchange interaction Hamiltonian with inhomogeneous exchange coupling 
strengths (Eq. (7)) and show that the entanglement gap temperature, Te, has a non- 
monotonic dependence on the exchange coupling inhomogeneity parameter a. The critical 
entanglement temperature, T^, obtained by using the susceptibility as an EW, has a mono- 
tonic dependence on a. We next determine how the entanglement gap temperature, T^, 
varies as a function of S in the cases of small spin clusters like a dimer, a trimer and a 
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tetrahedron. While Te increases with S in each case, the scaled entanglement gap temper- 
ature Ie decreases as S increases and goes to zero S ^ oo. The physical interpretation is 
that the entangled states have a small contribution to the total energy range in the limit 
of large S. The general applicability of this result for spin clusters without "all-to-all" spin 
couplings should be investigated. Toth |TT| has considered a Hamiltonian with "all-to-all" 
couplings between spin-| particles. The entanglement gap temperature Te is found to 
increase as increases but tE tends to a constant value as N becomes large. Bowling 
et al. have given examples of Hamiltonians describing bipartite systems for which tE 
increases without bound as the dimension of the Hilbert space associated with the subsys- 
tems increases. In our case, with increasing S', the Hilbert space of the system is enlarged 
but tE decreases as a function of S and goes to zero in the limit S oo. This is so since 
Te has a linear variation with S (Fig. 6) and Etot varies as S"^ in the large S limit. 

Lastly, we study a spin-1 dimer compound as an illustration of the quantum comple- 
mentarity relation. In experiments, the compound exhibits low-temperature magnetization 
plateaus. Our theoretical calculations reproduce these plateaus and further show that if the 
system is entangled, the magnetization plateaus coexist with the entanglement plateaus. 
Successive plateaus are connected by sharp changes in the magnetization and the amount 
of entanglement. The increase in one quantity is compensated by a decrease in the other 
quantity so that the complementarity relation is not violated. A large number of AFM 
compounds exhibit the phenomenon of magnetization plateaus [2I|. If these systems are 
entangled at the temperatures for which magnetization plateaus are observed, one can 
predict the coexistence of magnetization and entanglement plateaus in such systems. The 
Oshikawa, Yamanaka, Affleck (OYA) [28] theorem provides the condition for the occur- 
rence of magnetization plateaus in quasi-ld AFM systems. Magnetization plateaus have 
also been observed in a two-dimensional S=| AFM system SrCu2 {303)2, ^^^^ extending 
the scope for the applicability of the OYA theorem. It will be of interest to establish a 
connection between the OYA theorem and the quantum complementarity relation so that 
the conditions for the simultaneous appearance of the magnetization and the entanglement 
plateaus are clearly identified. 
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Appendix A: Eigenvalues and eigenvectors of linear chain 
tetramer 

The Hamiltonian describing the hnear chain tetramer is given in Eq. (7). The total spin 
of the tetramer is S"*"*. The first index in the subscript of an eigenvector refers to the 
eigenvalue and the second to 5**°*, the z-component of the total spin. 
gtot ^ 2 : 

£i = (I + f ) J (^1) 

^i,2 = |TTTT) 
^1,1 = HlittT) + inrt) + ITTiT) + ITTTi)) 

i^i,o = ^ (iTTU) + inti) + iTiiT) + urn) + unt) + liitt)) (^2) 

^1,-1 = I (ITiii) + litii) + liiTi) + liiiT)) 

V'l -2 = liiii) 

gtot ^ ^ . 

E2 = {-l + j)J {A3) 

V'2,i = i(iiTTT)-iTirT)-irnT) + iTTn)) 

^2,-1 = I (liTTT) - ITiTT) - ITTiT) + ITTTi)) (^4) 
V'2,o = ^(|TUT)-Um)) 

^3,1 = (liTTT) - ITTti) - (« - VT+a^) (ITiTT) - ITTiT)) 
^3,-1 = 1^ (ITiii) - liiiT) - (« - VTT^) (liTii) - liiTi)) {A6) 
i^3,o = i (iTTii) - liiTT) - - (ITiii) - liTiT)) 

E, = (-^ - ^^TT^) {A7) 

^4,1 = 4 (liTTT) - ITTTi) - {a + VT+ ^) (ITiTT) - ITTiT)) 
V^4,-i = i (ITiii) - liiiT) - (« + VTT^) (liTii) - liiTi)) (^8) 
^4,0 = it (ITTii) - liiTT) - (^ + VTT^) (iTiTi) - liTiT)) 

gtot ^ Q . 
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V'5,0 = ^ (ai ITTii) + h \nn) + ci Itiit) + litti) + ei litit) + /i liiTT)) (^10) 

-'V5 



V's.o = ^ (a2 ITTii) + b2 ITiTi) + C2 ITiiT) + ^2 liTTi) + e2 liTiT) + /2 liiTT)) (^12) 



where a,^h^a,^h^ 1, 61 = -f + 2^i " 2^ + ^ = ^i, &2 = -| - - ^ + ^ = 

^2, ci = -1 + I - 2^i - ^ + ^ = rf,, C2 = -1 + I + 2^i-^ + ^ = d, 
Ni, N2, N3, N4, N5 and Nq are the appropriate normalization constants. 

Appendix B: Eigenvalues and eigenvectors of the spin-1 
dimer 

The dimer Hamiltonian is given by Eq. (21). The basis functions are represented 
as with = ±1,0 and 5*1 = ±1,0. The eigenstates and the eigenvalues are 

described by ^„ m and A„ ^ where 'n' and 'm' refer to the total spin S'*"* of the dimer and 
its z-component respectively. 
5'*°* = 2 : 

02,±2 = |±1,±1) 

A2,±2 = 5-/ ± 2d ± 2S ^ ' 

02,±l = ;^(|±l,O)±|O,±l)) 

A2,±i = J + d±B ^ ' 



0: 



'2,0 



2 (^m 1 1 5 " 
^2,0 



■) + |- 

-f ±«! 



l,l)) + v^Ap|0,0) ^^3^ 
R 



not 



1 : 



ii,±i = ^(|±l,0)-|0,±l)) 
Xi±i = - J + d±B 



m 
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,,o = j-,{\l,-l)-\-l,l)) 
A2 ±1 = —SJ + 2d 



gtot ^ Q . 



00,0 = \ {Ap (|1, -1) + 1-1, 1)) - y/^Am |0, 0)) 



Ao,o — —^ + d — R 



where 



R 



{B7) 



R 
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